Diffusion-controlled annihilation A + B — > 0: The growth of an A particle island from a 

localized A-source in the B particle sea 
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We present the growth dynamics of an island of particles A injected from a localized A-source into 
the sea of particles B and dying in the course of diffusion-controlled annihilation A + B — > 0. We 
show that in the Id case the island unlimitedly grows at any source strength A, and the dynamics of 
its growth does not depend asymptotically on the diffusivity of B particles. In the 3d case the island 
grows only at A > A c , achieving asymptotically a stationary state (static island). In the marginal 
2d case the island unlimitedly grows at any A but at A < A* the time of its formation becomes 
exponentially large. For all the cases the numbers of surviving and dying A particles are calculated, 
and the scaling of the reaction zone is derived. 
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For the last two decades the reaction-diffusion system 
A + B — > 0, where unlike species A and B diffuse and ir- 
reversibly react in a d dimensional medium, has acquired 
the status of one of the most popular objects. This at- 
tractively simple system, depending on the initial condi- 
tions, displays a rich variety of phenomena and, depend- 
ing on the interpretation of A and B (chemical reagents, 
quasiparticles, topological defects, etc.), it provides a 
model for a broad spectrum of problems [1]. A crucial 
feature of many such problems is the dynamical reaction 
front - a localized reaction zone which propagates be- 
tween domains of unlike species. The simplest model of 
a reaction front, introduced more than a decade ago by 
Gain and Racz [2], is a quasi-one-dimensional model for 
two initially separated reactants uniformly distributed on 
the left side (x < 0) and on the right side (x > 0) of the 
initial boundary. Taking the reaction rate in the mean- 
field form R(x,t) — ka(x,t)b(x,t)(k being the reaction 
constant) Galfi and Racz discovered that in the long- 
time limit kt — > oo the reaction profile R(x, t) acquires 
the universal scaling form 



R = RfQ 



x f 



(1) 



where xj oc t 1 ! 2 denotes the position of the reaction zone 
center, Rf oc t~P is the height, and w oc t a is the width 
of the reaction zone. Since then, much work has been de- 
voted to studying this problem by different approaches. 
It has been shown [3 - 7] that the mean-field approxi- 
mation can be adopted at d > d c — 2 (with logarithmic 
corrections in the 2d case), whereas in Id systems fluc- 
tuations play the dominant role and the explicit form of 
R remains unknown. Nevertheless, the scaling law (1) 
takes place in all dimensions with a = 1/6 at d > d c and 
a = 1/4 at d = 1, so that at any d the system demon- 
strates a remarkable property: on the diffusion length 
scale Lb oc t 1 / 2 the width of the reaction front asymptoti- 
cally unlimitedly contracts: w/ Lb — > as t — > oo . Based 
on this property a general concept of the front dynam- 
ics, the quasistatic approximation (QSA), has been de- 
veloped [3, 5, 8, 9]. The QSA consists in the assumption 



that for sufficiently long times the dynamics of the front 
is governed by two characteristic time scales. One time 
scale, tj = — (dlnJ/dt) -1 , controls the rate of change in 
the diffusive current J = J a = Jb of particles arriving 
at the reaction zone. The second time scale, tf oc w 2 /D, 
is the equilibration time of the reaction front. Assum- 
ing tf/tj <C 1 from the QSA in the mean-field case with 
equal species diffusivities Da,b = D it follows [3,8] 



R 



f 



J/w, w - (D 2 /Jk) 1 / 3 , 



(2) 



whereas in the Id case w acquires the fc-independent form 
w ~ (D/ J) 1 / 2 [3,5,9]. The most important feature of the 
QSA is that w and Rf depend on t only through the 
boundary current J(t), which can be calculated analyt- 
ically without knowing the concrete form of Q. On the 
basis of the QSA the general description of the system 
A + B — > with initially separated reactants have been 
obtained for arbitrary diffusion coefficients [10]. These 
results are in full agreement with extended numerical cal- 
culations and have been generalized recently to the case 
of nonmonotonic front motion [11]. 

The purpose of present paper is to apply the QSA to 
the long-standing problem of growth of an A-particle is- 
land from a localized A-source in the uniform B particle 
sea. This important problem was first analyzed by Lar- 
ralde et al. [12] for the special case of static sea (Dg = 0). 
Assuming that diffusing A particles are injected at a sin- 
gle point into a reactive d-dimensional substrate B and 
instantaneously react with B upon contact, Larralde et. 
al. have studied the growth dynamics of the reacted re- 
gion radius, r/(f), and the number of dying and that 
of surviving A particles. Considering the reaction front 
dynamics as a Stefan problem, they have, in particu- 
lar, shown that at any source strength r / asymptotically 
grows by the laws oc (tint) 1 / 2 , oc t 1 / 2 and oc t 1 / 3 at 
d =1,2 and 3, respectively. Subsequently, those results 
were generalized to the cases of imperfect reaction [13] 
and diffusion with a bias [14], however, as in [12], the B 
particles were always presumed "frozen". In this Rapid 
Communication we present a theory of growth of a d- 
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dimensional A-island for the physically most important 
situation when both A and B particles are mobile. In 
the framework of the QSA we first consider the simplest 
"standard" case with equal species diffusivities, and then 
we extend the obtained results to the case of arbitrary 
nonzero diffusivities, thus revealing a rich general picture 
of the island growth for d=l, 2, 3. 

Let at t > particles A be injected with a rate A at the 
point r = of the uniform c?-dimcnsional sea of particles 
B, distributed with a density p. Particles A and B diffuse 
with nonzero diffusion constants D AB an d upon contact 
with some nonzero probability annihilate, A + B — > 0. In 
the continuum version such process can be described by 
the reaction-diffusion equations 

^=D A V 2 a-R + A6(r), ^ = D B V 2 b-R (3) 

with the initial conditions a(r, 0) = 0, b(r, 0) = p, and the 
boundary condition b(r — > oo,t) = p. Here a(r,t) and 
b(r, t) are the mean local concentrations of A's and B's 
which, by symmetry, we assume to be dependent only on 
the radius, and R(r, t) is the macroscopic reaction rate. 
Equal species diffusivities. — To simplify the problem es- 
sentially we will first assume, as usually, D A = D B = D. 
The initial density of the sea, p, defines a natural scale 
of concentrations and a characteristic length scale of the 
problem - average interparticle distance I = p~ x l d . So, 
by measuring the length, time and concentration in units 
of £,£ 2 /D and p, respectively, we introduce the dimen- 
sionless source strength A = A£ 2 /D and the dimension- 
less reaction constant k = k / '^( d ~ 2 ) D . Defining then the 
difference concentration s(r, t) = a(r, t) — b(r, t) we come 
from (3) to the simple diffusion equation with source 

ds/dt = V 2 s + A(5(r) (4) 

at the initial and boundary conditions s(r, 0) = s(r — > 
oo, t) = —1. According to (4) in the course of injec- 
tion in the vicinity of the source there arises a region 
of A particle excess, s > 0, which expands with time. 
The central idea of the paper is that, by analogy with 
the Galfy-Racz problem, a narrow reaction front has to 
form at this region boundary, of which the law of mo- 
tion, rf(t), according to QSA, can be derived from the 
remarkably simple condition s(rf,t) = 0. Then, under 
the assumption that on the scale rj the front width w 
can be neglected w/tf <C 1, i.e., setting that a = s, b = 
at r < 77 whereas a = 0, b = \s\ at r > r/, the number 
of surviving N A {t) and that of dying N x (t) A particles 
is immediately derived from the condition 

N A = At - N x = fi d / f s{r, ty^dr (5) 
Jo 

with fii = 2, fi 2 = 27T , fi3 = 47r. By calculating in the 
limit wjrj <C 1 the current of particles in the vicinity 
of the front, J = — ds/dr\ r=rf , from Eqs.(2) one can 
easily obtain the law w(t) and define, in the end, a self- 
consistent condition of crossover to a quasistatic scaling 



regime (1). We start with an analysis of the behavior of 
rf(t),N A (t) and N x (t) for each dimension separately. 
Id case. — In Id the solution to Eq.(4) has the form 

s(r, t) = {V>H)ieric{r/2y/t) - 1, (6) 

whence, according to condition s(r/, t) = 0, the equation 
of motion of the reaction front center, r/(i), is 

ierfc(r / /2Vt) = l/VA 2 ^, (7) 

where ierfc(C) = erfc(v)dv = e - ^ /^/n — £erfc(£). 
From (6), (7) it formally follows that an excess of A par- 
ticles forms in a time t c = ir/\ 2 . It is, however, clear 
that a continual approximation comes into play at times 
t 3> max(l,l/A), therefore at early island formation 
stages one can distinguish two qualitatively different is- 
land growth regimes: i) A <C 1, when the island formation 
proceeds under conditions of death of the majority of in- 
jected particles and ii) A 3> 1, when the island forms 
long before the beginning of intensive annihilation. Let 
us consider first the limit A <C 1. In this limit an interval 
between injection acts, 6t\ = 1/A, is quite large, there- 
fore in the case of perfect reaction each injected particle 
dies long before the next one appears until the distance to 
the nearest sea particle, ~ At — \ft, becomes comparable 
with the characteristic diffusion length y/8t\. Whence for 
the time of beginning of the injected particles accumula- 
tion we find tb oc t c oc A~ 2 , which reveals the sense of t c . 
Assuming e = (t — t c )/t c < 1 we have Q = r//2\/t <C 1 
and hence from (5)- (7) we find 

r f ~ e/A, N A ~ e 2 /2X, N x ~ tt(1 + e)/A, 

whence it follows N A /N X cx e 2 . Calculating from (6), (7) 
the current J ~ A(l — e\/ir)/2 and assuming k not 
to be too small [9](k > VA) from the fluctuation law 
w ~ J -1 / 2 ~ l/\/A we find w/rj ~ \f\je ~ \j\[W A 
and tf/tj cx A 2 . Thus, the condition of crossover to the 
regime of quasistatic front is e > \/X(N A 3> 1). Defin- 
ing a minimal island by condition w/rj ~ N A ~ 1 for 
the island formation time we have ~ \/A and hence 
tb ~ t c (l + \/A) and rb ~ 1/VA- In the long-time limit 
T = t/t c 3> 1 from (7) it follows Q » 1 and we can 
rewrite (7) in the form 2e^fQ- = \pf whence we obtain 
exact asymptotics 

r f = V2~tWT{l -In lnT/lnT +•••), (8) 

and from (5), (6), and (8) we find N A — \t[\ — 
0(y/\nT/T)],N x = r/(2 + (J 2 + •••) = V%tWT. Con- 
sider now the limit A » 1. It is evident that in this 
limit a multiparticlc "cloud" forms long before the be- 
ginning of noticeable annihilation, therefore the stage of 
developed reaction (8), t 3> 1(3> t c ) is preceded here 
by the stage of purely diffusive expansion of the cloud, 
1/A <C t <C 1. The statistical theory of diffusive cloud 
expansion has been developed recently in the work [15]. 
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According to [15] in the "collective" regime (At 3> 1) 
the "radius" of an Id cloud grows (in our units) by the 
law r + = y/At ln(Ai). Remarkably, r + and r/ "join" 
exactly in the front formation zone (t ~ <~ r/) 

whereas at t > lr/ begins to grow slower than r+, as 
it has to. So, forming in qualitatively different regimes 
from p = N X /N A > 1(A « 1) to p « 1(A > 1) the 
Id island at any A crosses over to the universal growth 
regime (8) with an unlimited decay of the dying parti- 
cles fraction p oc ^/ln T/T — > 0. It remains for us to 
reveal conditions of quasistaticity of the front (8). Ac- 
cording to (6) in the limit i,T>l current J ~ \J\nT~/i. 
Thus assuming k not to be too small (k > we find 
w ~ J-V2 _ (i/lnT) 1 / 4 whence w/r f ~ (Hn 3 T)-V 4 
and — (ilnT)~ 1/2 . As C/ > 1, the conditions 

w/rf,tf/tj <C 1 ought to be supplemented by a more 
strict requirement of equality of currents at both front 
sides w < C = \ r =r f = r f /20. Calcu- 

lating w/C <~ (InT/t) 1 / 4 we arrive at the requirement 
i 3> max(t c , InT). 

2d case. — In 2d the solution to Eq.(4) has the form 

s(r, t) = -(A/47r)Ei(-r 2 /4i) - 1, (9) 

whence, according to condition s(rf,t) = 0, the equation 
of motion of the reaction front center, r/(t), is 

Ei{-r 2 f /tt) = -4tt/A. (10) 

Here Ei(— Q = — dve~ v /v is the exponential integral 
which has the asymptotics Ei(— £) = \n(-f() + - ■ ■ at ( <C 1 
(7 = 1.781...) and -e~ c /C, + ■ ■ ■ at C > 1. From (10) it 
follows that rf grows by the law 

r f = 2\/cd, (11) 

where a is the root of the equation Ei(— a) = 
— A*/A, A* = 47r and has the asymptotics a = e~ x '/ x /j 
at A <C A* and a = ln(A/A*a) at A > A*. From (5), (9), 
(11) it follows 

N A = \t(l-e- a ), N x =Xte- a . (12) 

We conclude that in 2d the island growth rate a and the 
dying-to-surviving particles ratio p do not vary in time: 
at large A ^> A* the majority of particles survive, p ~ 
In A/ A, whereas at small A <C A* the majority of particles 
die, p ~ e A */ A . The most interesting consequence of 
(10) consists in the exponentially strong decrease of the 
growth rate in the region A < 1. Defining a minimal 
island through the condition Na <~ 1 for its formation 
time at A < 1 we have % ~ e A */ A /A whence it is seen 
that at A < A» the island growth is actually suppresed. 
Calculating current J = X/X*e a \/ai for the scaling of 
the reaction zone from (2) we find w ~ (t/tw) 1 ^ ,t w — 
(nX/e a y/a) 2 . At A C A» this yields yftf/tj < w/rf ~ 
(tb/nt) 1 / 3 whereas at A 3> A* we have \ftfjtj Cra/£~ 
(In X/nt) 1 / 3 . Thus, crossover to the quasistatic regime 



occurs at times nt 3> % and nt ^> In A, respectively (note 
that k < 1, being ~ 1 for perfect reaction). 

3d case. — In 3d the solution to Eq.(4) has the form 

s(r, t) = (A/47rr)erfc(r/2Vt) - 1, (13) 

whence, according to condition s(r/, t) = 0, the equation 
of motion of the reaction front center, rf(t), is 

erfc(r//2Vi) = 47rr / /A. (14) 

From (14) it follows that Q(t) decreases indefinitely 
so that at large t/t s » 1 the front radius, by the 
law rf — r s [l — 0{y/t s /t)] with the characteristic time 
t s = (A/A*) 2 , reaches a stationary value 

rf(t/t s — » 00) = r s — A/A*. (15) 

According to (5), (13), (14) in this limit N A = ^-r 3 [l - 
0{yfQi)\,N x = Xt[l - 0{t s /t)] therefore, in contrast 
to the Id case, at any A all the injected particles die. 
The steady-state current J s = A* /A, whence, accord- 
ing to (2) w s ~ (A/A*k) 1/3 and w s /r s ~ (A^/A^/k) 2 / 3 . 
Defining a minimal stationary island through the con- 
dition w s /r s ~ 1 we conclude that in the 3d case the 
island forms only when the injection rate exceeds a criti- 
cal value A c ~ X*/y/n. A maximal value of k, attainable 
in the perfect reaction limit, is k p ~ a/l {a being the 
size of particles), therefore A c > A* 3> 1. By rewriting 
Eq.(14) in the form erfc(£/)/C/ = 2y/t/t s one can easily 
see that at high injection rates A/ A* 3> 1 the station- 
ary stage (t » t s ) is preceded by an intermediate stage 
1 <C t <C t s wherein the island grows by the law 

r f = y/2t\n.{t./t){\ - ln(vMA> + •••), (16) 

where w = \n(t s /t). According to (5), (13), (16) at 
this stage N A = Xt[l - O{uo 3 ' 2 y/i/U)],N x = ^-rj[l + 
0(w -1 )] and, therefore the majority of particles have yet 
survived, p <~ oj 3 / 2 \ftjt s <C 1. Calculating the current 
J <~ y/uj/t we find from (2) w ~ (t/n 2 ^) 1 / 6 whence 
y/tf/tj < w/C - (uj/nt) 1 / 3 . Thus, the formed front 
condition reads t ^> \n(t s /t)/ n. According to [15] the ra- 
dius of a 3d cloud, which expands in t he absence of reac- 
tion, in our units has the form r + = \J 2t\n[(Xa / l)H / Air]. 
Comparing r+ and r/ suggests that r/ begins to lag be- 
hind r+ at times t 3> t/a ~ 1/k p in remarkable agree- 
ment with the above estimation. 

To sum up the above, as key results we distinguish 
a,) self-consistent analytic description of the island growth 
since the moment of its formation and b) anomalously 
slow island growth at A < A* in 2d and complete suppres- 
sion of its growth at A < A c in 3d, which contrast sharply 
with growth asymptotics in the static sea [12]. 

Arbitrary ratio of diffusivities. — Let us now extend the 
analysis to a general case < V = Db/Da < 00, nondi- 
mensionality with respect to D = Da being retained. We 
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present final results here, a detailed discussion of which 
will be given elsewhere. 

[d=l]. Comparing our results with the results of 
Ref.[12] we find that long-i asymptotics (LTA)(8) for 
T> = 1 converges to the LTA for I? = 0. We thus con- 
clude that the growth of the Id island does not depend 
asymptotically on the diffusivity of B particles. This 
conclusion is the consequence of the evident fact that at 
rf/y/Vt 3> 1 particles B can be regarded as effectively 
static, so as rf/\/i oc Vint — > oo the LTA of the is- 
land growth at any < V < oo must converge to the 
LTA for static sea. In the interval < T> < 1 the time 
of crossover to the LTA does not alter appreciably al- 
though the time of A particle accumulation, t h , shifts 
considerably at small A (comparing Xt — \pDt with y/8t\ 
we find that in the interval V <C VX <C 1 the value of 
tb oc A~ 3 / 2 does not depend on V whereas at \/A <C V 
it grows with V by the law % oc 2?/A 2 ). At T> 3> 1 the 
time of crossover to the LTA becomes exponentially large, 
t ^> eP /A 2 , therefore the transient dynamics in this limit 
requires to be specially considered (note that for this case 
at small \JXjD <C 1 we again find % oc V/X 2 ). 

[d=2] . In the 2d case the solution of the problem with a 
source possesses a remarkable property: s(r,t) — f(-^j)- 
Using this property and assuming w/r/(t —> oo) — > it 
is easy to check that the solution of Eqs.(3) has to read 
a = -(^)Ei(-C 2 ) - A, b = at r < r f and a = 0, b = 
1 + BEi(-( 2 /V) at r > r f with C = -^(A,B = const). 

Equalizing A's and B's currents at r = rf we come to the 
laws (11), (12) with the exact equation for a at arbitrary 
V: 

Ei(-a/X>) = -{AirV/X)e a{v - 1)/T) . (17) 
In the limit V < 1 from (17) it follows 

C Ve- X * v ' x h, A/A, « V 
a ~ I A/A*, V < A/A* < 1 

[ ln(A/A*a), A/A* > 1. 

In the opposite limit V >• 1 from (17) it follows 

< x> e -A.A/ 7) A/ A* < In -1 © 

ln[Aln(X>/a)/A*], In" 1 V < A/A* < e p 
I ln(A/A*a), A/A* > e p , 

where A = A/2?. Thus, at strong difference of dif- 
fusivities we find three characteristic growth regimes: 
a <C a- = min(l, 2?)(I), a_ <C a <C a + (ll) and 
a 3> a+ = max(l,2?)(III). In regime I the majority of 
particles die, and the island growth does not depend on 
Da ■ In regime III the majority of particles survive, and 
the island growth does not depend on Db- In interme- 
diate regime II at V <C 1 the majority of particles die 
and the island growth does not depend on Db whereas 
at V 1 the majority of particles survive and the island 
growth depends on the both diffusivities. Essentially, 



that at any finite T>, as A decreases, the island growth 
crosses over to the regime of " exponential suppression" I 
which disappears only in the limit T> — > 0. 

[d=3]. One can easily check that in the general case 
T> 7^ 0, like in the case T> = 1, the 3d island has to 
reach asymptotically a stationary state. Indeed, assum- 
ing ^ — > in the steady state limit t — > oo from Eqs.(3) 
we find a s = 6(r s - r)V(?f - 1), b s = 0(r - r a )(l - ^), 
where 0(x) is the Heaviside step function and stationary 
radius 

r s = A/A*2? = M 2 /4irD B 

does not depend on Da- For the particle number we find 
NX = ^fVrl ~ (A/A*) 3 /X> 2 whence we obtain the lower 
island formation boundary A/ A* 3> 2? 2//3 and conclude 
that a mean island density < a > s = D/2 does not depend 
on A: the island is always concentrated (with respect to 
the sea) at V 3> 1 and is always rarefied at V <C 1. 
According to [10] for V ^ 1 w s ~ (£j) 1/3 ~ (if) 1/3 
whence it follows w s /r s ~ (r 2 /t) -1 / 3 . Thus the necessary 
condition for the 3d island formation reads 

A > A c - (A*2?/V^)max(l, \j k/V 2 / 3 ). 

We have been unable to describe analytically the com- 
plete kinetics of crossover to the steady state for arbi- 
trary T> ^ 1. However, in the interval < T> < 1 in- 
termediate asymptotics of the 3d island growth can be 
revealed based on simple arguments. Indeed, at V < 1 
in the limit \Tb~i 3> 77 we have 77 <~ r s . In the oppo- 
site limit y/vi «rj the sea is effectively static, therefore 
the island ought to grow by the law r/ ~ (At) 1//3 [12]. 
From the both conditions for the crossover time it fol- 
lows t s oc X 2 I'D 3 so that at t < t s the particle number 
grows by the law Na ~ (A 5 / 2 i) 2 / 3 , and the density de- 
cays to < a > s by the law < a >^ (A 2 //;) 1 ' 3 . Clear that 
at A ^> 1 this asymptotics is preceded by the law (16) 
(crossover time U <~ A 2 ) whereas at A c <C A <C 1 the 
island formation time tb oc A -5 / 2 . In the limit 2? — > it 
follows A c — > 0, r s — ► 00, t s —* 00 and we are coming back 
to the unlimited island growth at arbitrary finite A. 
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